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' We develop a simple computational model for cell boundary evolution in plastic deformation. 

We study the cell boundary size distribution and cell boundary misorientation distribution that 
experimentally have been found to have scaling forms that are largely material independent. The 
' cell division acts as a source term in the misorientation distribution which significantly alters the 

scaling form, giving it a linear slope at small misorientation angles as observed in the experiments. 
We compare the results of our simulation to two closely related exactly solvable models which exhibit 
scaling behavior at late times: (i) fragmentation theory and (ii) a random walk in rotation space 
CZ2 , with a source term. We find that the scaling exponents in our simulation agree with those of the 

theories, and that the scaling collapses obey the same equations, but that the shape of the scaling 
functions depend upon the methods used to measure sizes and to weight averages and histograms. 
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INTRODUCTION 

a 

After significant plastic deformation, the dislocation tangles in crystals often organize themselves into sharp walls 
separating nearly dislocation-free cells; the crystallographic axes rotate slightly across each cell boundary. These cells 
undergo refinement (become smaller) under increased deformation, and recent experiments^, El indicate that both 
the misorientation angles and the cell sizes have power-law scaling with material-independent scaling forms for the 
probability distributions. We introduce here a simple model of cell division and rotational diffusion which exhibits this 
type of scaling, and which provides insights into the origins for the experimental scaling distributions. In particular, 
1^ ' we argue that cell division (driving the refinement) is responsible for the linear growth of the misorientation scaling 
, distribution at small angles. 

Cell boundaries are distinct from grain boundaries in that their misorientation angle across them is small (at most 
a few degrees) and they form in a nonequilibrium process, typically at temperatures where diffusion is not relevant 
(so, for example, the impurity segregation characteristic of many grain boundaries is not observed at cell boundaries). 
As deformation proceeds, the cell structure refines (the average cell size L av becomes smaller), and the average cell 

S misorientation angle av grows. 
The cell boundaries are separated into into two classes. An early workQ called them "ordinary cell walls" and 
"dense dislocation walls"; later authors have called them GNBs ("geometrically necessary boundaries") and IDBs 
("incidental dislocation boundaries"). The GNB's typically align roughly parallel to one another, have larger misori- 
entation angles, and are longer, often forming the boundaries of two or more cells. 

(The term "geometrically necessary" is unfortunate. Geometrically necessary dislocations are those required to 
mediate macroscopic strain gradients and rotation gradients, as distinguished from geometrically unnecessary dis- 
locations whose Burgers vectors cancel out on long length scales. All cell walls are associated with small relative 
rotations between cells, and are hence composed of geometrically necessary dislocations on scales comparable to the 
cell sizes. On longer scales, far from building up large macroscopic rotations, the rotations mediated by neighboring 
geometrically necessary boundaries tend to cancel 13J, leading to little or no long-range rotation gradient. Hence, as 
cell boundaries, the GNB's are most akin to the geometrically unnecessary dislocations.) 

Hughes et al.fy\^ studied the distribution functions for these two types of cell walls, and found a simple scaling 
behavior, largely independent of material. In particular, for the GNBs, 9 av ~ e 2 / 3 and L av ~ e~ 2 / 3 , while for the IDBs 
Sav ~ e 1 / 2 and L av ~ e -1 ^ 2 , where e is the magnitude of the net plastic strain. Moreover, data for several materials 
and different strain amplitudes all collapse onto apparently universal scaling curves p m i s and p s i ze when rescaled to 
the average angle: 

p(O) = 9- v 1 p mis (0/9av) (1) 
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and 



P( L ) = Ljp sizc (L/L av ). 
We will study the scaling of these probability distributions p(9) and p{L) using a simple model. 



(2) 



MODEL 



How much of this apparently universal behavior can be captured in a simple model of cell rotation and refinement? 
The model we propose is one in which cells become smaller by subdivision (leading to a fragmentation theory for the 
size distribution), and undergo random angular reorientations as the strain increases. Our model does not incorporate 
the anisotropy in the external strain field, and so has nothing to say about how cell structure morphology might change, 
say, between tensile and rolling deformation or as the crystalline orientation changes. One should view our model as 
a caricature of the real system; our results suggest that the experimentally observed scaling behavior may be generic 
to any microscopic mechanism which fragments and randomly reorients cells. 

Our computational model starts with one large cubical cell. We assume a cubic crystal, with initial crystalline axes 
aligned with the axes of our cube, so the initial orientation is described by a rotation matrix i?(0) equal to an identity 
matrix. The dynamics of our model incorporates two pieces: rotational diffusion and cell splitting events. 
Rotational Diffusion. The orientation of each cell a undergoes a simple random walk in rotation angle space, with 
strain increments playing the role of the time step. It is convenient to write the current orientation R a {t) = exp(n-J), 
where the matrix = e^fc (with e^-fe the totally antisymmetric tensor) generates an infinitesimal rotation about the 
j th axis. Since the cell boundary misorientation angles experimentally are small (around a degree or so), we may 
expand the exponential in this expression, 



(Large angle corrections are discussed in |4| and |{J.) In this approximation, diffusion in the manifold of crystalline 
orientations can be written as an ordinary diffusion equation in the three-dimensional coordinate n. If we assume e(t) 
is a monotonically increasing strain, then the three-dimensional probability distribution of grain orientations A(n) 
evolves according to the equation 



where D is the "orientational diffusion constant" and the Laplacian V 2 = d 2 /dn 2 + d 2 /dn 2 + d 2 /dn 2 . The random 
walk described by this diffusion equation is implemented numerically by adding a Gaussian random vector to n with 
components of root-mean-square length V2~DAe whenever the strain for the cell is incremented by Ae. 

Cell Splitting Events. Our model, for simplicity, divides cells only along planes perpendicular to one of the 
crystalline axes. Thus our cell structure is composed of rectangular parallelepipeds. The rate of cell division in our 
model depends only on the current size and shape of the cell, and not on its environment. There are several different 
physical mechanisms that might be responsible for cell division. Broadly speaking, we classify them by dimensional 
analysis: there are mechanisms that divide cells at a rate proportional to their current volume V, their current surface 
area S, their current perimeter P, or at a uniform rate independent of the current size {7 = 1. Once a cell has been 
chosen to split, we must choose an axis and a position along that axis to place the new cell wall. To keep our aspect 
ratios reasonable, we've chosen the probability of splitting along a given axis proportional to the length of the cell 
along that axis. The position of the new cell wall along the split axis is chosen randomly in all cases. The two cells 
formed by splitting inherit their parent's orientation: the new cell walls thus start out at zero misorientation angle, 
which will be important when we study the misorientation angle distribution. 

Our model for cell splitting is closely related to a well-studied model of fragmentation^, In fragmentation 
theory, the splitting rate is assumed a function of the volume, so for example a cell of volume V could fracture with 
a rate AV 1 . In our problem, 7 = and 7=1 correspond precisely to uniform and volume cell splitting rates, while 
7 = 1/3 and 7 = 2/3 approximately correspond to perimeter and area splitting. The evolution law in fragmentation 
theory corresponding to our model is easily seen to be|| 



R = 



exp(n • J) 




(3) 



dA(n)/de = -L»V 2 A(n), 



(4) 




(5) 
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where c is the concentration of fragments with volume V . To relate this to experimental measurements (equation [3J, 
which produce probability distributions of lengths rather than concentrations of volumes, we can change variables 
from V to L — V 1 / 3 . The probability distribution of lengths with all cells weighted uniformly is 




FIG. 1: Simulated Cell Morphologies. Diagonal cross section of the cell morphologies in two model simulations, in a central 
plane perpendicular to the 111 axis. LEFT: Area splitting rate. RIGHT: Uniform splitting rate. The area splitting yields 
rather uniform cell sizes, while the uniform splitting yields an enormous range of cell sizes and a fractal morphology. The colors 
are chosen to represent the rotations of the crystalline axes of the individual grains. The original orientation n = is colored 
gray; n x , n y , and n z are mapped respectively onto deviations in red, green, and blue with a scale-factor chosen to saturate at 
the largest rotations. 

The histograms produced by our simulations are a result of three choices. First, there is a choice in how we define 
the size (length L) of the cell. In computing the averages and histograms from the simulation data, we typically 
define the size of a cell to be its length along any one of the axes: all three lengths are incorporated into the averages 
and histograms. This definition of size corresponds to that used in experiment. Alternatively, in order to compare 
the histograms from the simulation to fragmentation theory (which keeps track only of the volumes, not the shapes, 
of the cells), we can define the size of the cell as the cube root of the volume. Second, there is a choice in the 
splitting dynamics as discussed above: cells can divide at rates proportional to their volume, surface area, perimeter, 
or at a uniform rate. In fragmentation theory, the rate of splitting is proportional to V 1 as discussed below: thus 
fragmentation theory is exact for our simulation with volume splitting (7=1) and uniform splitting rates (7 = 0), 
but does not directly apply to the perimeter and surface simulations, whose splitting rates depend upon the shapes 
of the cells as well as their volumes. Finally, one must address how to weigh the contributions of different cells in 
the probability distribution. For example, an experiment which measures cell sizes by taking an XY cross section 
and then weighting each observed cell equally in the average is effectively weighting the three dimensional cells by 
their extent in the Z direction (roughly weighting each cell by its perimeter or by L). We compute the averages and 
histograms from the simulation data by weighting each cell uniformly or by its volume, surface area, or perimeter. 
We will soon see that the scaling exponents for the average size depend on only the splitting dynamics, 7, not on 
the measurement of size or the weight of the distribution. We will also see that the shape of the scaled probability 
distributions changes with different measurements of size or weights, but the distributions scale nonetheless. 

If we define the size of the cell as the cube root of the volume, simulations that split cells at a uniform rate or at a rate 
proportional to the volume of the cell produce histograms that agree well with those given by fragmentation theory 
(equation . However, simulations that split cells at rates proportional to area and perimeter produce histograms 
that are shifted from equation JJjJ . We focus here on the case where the rate is proportional to the surface area S of 
the cell (Area) and the case where the rate is independent of the size and shape of the cell (Uniform). 
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RESULTS 
Morphologies 




FIG. 2: Mean Cell Sizes versus Strain. We measure the mean cell sizes, with each cell contributing to the average an 
amount given by volume, area, perimeter, or uniformly independent of its size. As noted in the text, the perimeter scaling 
corresponds effectively to a common experimental procedure of taking histograms. (Don't confuse, say, area weighted mean 
sizes and area weighted splitting rates. The different curves on each graph are different measurements, the different graphs 
are different dynamics.) We also include the mean size with the length measured as the cube root of the volume, with cells 
contributing uniformly, independent of their sizes. LEFT: Area splitting. The various measures are equivalent up to an overall 
scale factor for the model where the cell division rate is proportional to the surface area of the cells, and the cell size varies as 
This equivalence is because the distribution of sizes is peaked about a typical size scales (figure^, and can be derived 
from the scaling form for the cell size distribution. RIGHT: Uniform splitting. The various measures are very different for the 
case of random splitting, since the cell structures are fractal with cells of all sizes. 

Figures show the cell morphologies from the two simulations. Area splitting shows a fairly uniform density of 
cell sizes: this is characteristic also of the other size-dependent cell division rates. Uniform cell splitting rates produce 
a broad range of sizes: most cells are very small, so most cell divisions subdivide very small cells. Indeed, as we 
shall discuss below, the uniform model is at a critical point in fragmentation theory (the "shattering" transition^). 
Experimentally, there does not seem to be a consensus on whether the cell size distribution is fractal|9j or whether it 
has a more traditional scaling distribution with a characteristic size which shrinks with timeQ- 



Cell Size Scaling 

Figures (J2J shows the mean cell size as it evolves with increasing strain in our model, under area splitting and 
random splitting. In the figure, we show a count average where each cell contributes equally, a perimeter average 
where cells are weighted by their perimeter, and area and volume average. Each of these has the size of the cell 
defined as the distance between parallel cell walls. The fourth average has the size defined as the cube root of the 
volume and is weighted uniformly. 

The mean size for the cells in the Area splitting model scales with e -1 / 2 . Assuming fragmentation theory (equation^ 
and scaling (equation [3 which can be derived from fragmentation theory^ for 7 > 0), it is easyJlTj to derive the 
power- law relation L av ~ e -1 / 37 . Thus the experimentally observed scaling of the IDB sizes suggests a cell splitting 
rate proportional to the cell surface area (left side of the figures), while the scaling L [ £ NB) ~ e~ 2 / 3 of the GNB 
dislocations would suggest a mechanism with a cell splitting rate scaling as the cell volume to the 1/2 power. (See, 
however, part (2) in the conclusion.) 

Figures © show the histograms of cell sizes at the end of our simulation, with the five weights discussed above. 
Naturally, for example, if cells are weighted by their volume there is more weight in the histogram at larger sizes. 

We first consider the left-hand panel, showing the simulation results with dynamics which split cells proportional 
to their surface area. The Area splitting histograms, taken at different e during our simulations, do indeed rapidly 
converge to those shown on the left panel of figure (hence validating the scaling form equation (J2J). However one 
defines the size L or how one weights the contribution of each cell, the histograms collapse onto scaling functions. 
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FIG. 3: Cell Size histogram. Histograms of cell sizes at the end of our simulation, rescaled to the average cell size (figure 
|5J according to equation Again, there are different histograms depending upon how the cells are weighted in the average. 
LEFT: Area splitting. The dashed curve is the fragmentation theory prediction|| for the closely related problem where cells 
split at a rate proportional to volume V 2//3 . RIGHT: Uniform splitting. Notice that most of the cells (Counted) are very small. 



Fragmentation theory|8| calculates the scaling function explicitly that describes the distribution of sizes at late 
times. In particular, for Area splitting (7 = 2/3), the scaling function in (J2J for a uniformly weighted distribution is 

p size (x) = 32x 2 exp(-4a; 2 /7r)/7r 2 (7) 

shown as the dashed line on the left in figure As expected, it does indeed agree well with the size distribution 
where L is defined as the cube root of the volume. The shift between these two curves is due to the difference in 
dynamics: our simulation splits cells at a rate proportional to the actual area of the cell, not by V 2 / 3, . For Volume 
splitting, with L defined as the cube root of the volume, the histograms from the simulation agree with those from 
fragmentation theory as expected. 

The histograms which measure the length, width, and height of each cell are broader than those that measure the 
cube root of the volume. Cells with large aspect ratios will contribute one or two dimensions that are smaller than 
the cube root of volume and one or two that are larger. Notice that in figure © the cell size probability density does 
not vanish at zero size L except for statistical weightings that involve the total volume. In our model, a subdivision 
occurs with equal probability at all thicknesses, so the probability density at zero thickness is finite. If the weights 
of the cells in the average is by the total volume, of course, then the thin cells contribute vanishing weight so the 
histogram goes to zero. 

Consider now the right-hand panel of figure (0 , with simulation results for dynamics which split cells at a uniform 
rate independent of their size. Notice first that the scales are logarithmic: there is an enormous range of cell sizes, 
with peak probabilities (by most measures) at very small sizes. The Uniform splitting model is at a critical point in 
the parameter 7 = 0, the shattering transition Q, beyond which (7 < 0) there an infinite dust of zero-size particles. 
At this critical point, there are exact solutions for the cell size distribution (dashed curve shown in figure ©, to be 
compared with the L — V x l z simulation curve). This exact solution does not have the scaling form of equation J5J. 
We have not been able to find a generalization of the scaling form which collapses the distribution at the critical 
point, but the scaling variable (typical size) must shrink exponentially as e grows. One hint for why finding the 
scaling function may be difficult is that the system appears not to be self-averaging: the relatively good agreement 
between the theory and simulation shown on the right of figure © is not as true with other random number seeds, 
with fluctuations of an order of magnitude away from the theory. (The average over many seeds does agree with the 
calculated form.) While there is no scaling form solution to equation JSJ for 7 = 0, we have found a family of formal, 
non-normalizable solutions of the form c(V : e) = e _ ( 1 + 2 /^) € I/' 3 . 



Misorientation Scaling 

To study the misorientation angles, we need a formula for the misorientation angle 9. Imagine a cell boundary as 
in figure QJ. The rotation R^R^ 1 takes the crystalline axes of right-hand cell to the orientation on the left-hand side 
of the boundary. As before, we can write this rotation as exp(n • J), and correspondingly write Rl and Rr in terms 



FIG. 4: Geometry of a Cell Boundary. 



of rii and n^. Because the rotation angles are small across cell boundaries, n ~ riL — nR. The misorientation angle 
is 

T^RlR^ 1 ) = 1 + 2cos6» ~ 3 - 9 2 = 3 - n 2 

6~ ^nf + nj + n 2 = |n| (8) 

Define S(n)d 3 n to be the area of cell boundary with misorientation matrix R = cxp(n- J). Since this distribution for 
our model is symmetric under rotations, equation (jHJ implies that the probability distribution for the misorientation 
angle 

p(0) = p{\n\) = WS(n)A4 tot (9) 

where A to t — J d 3 riE(n) is the total cell boundary area. Notice that the (one-dimensional) probability density for 
small angles 9 is 4tt9 2 /A tot times the (three-dimensional) cell boundary area at one of the rotations n corresponding 
to 9 = |n|. This is of course because the number of possible rotations grows with misorientation angle, like the 
area of a sphere in rotation space. This has the important consequence of making p{9) vanish at 9 = 0; this 
reflects not some special physics which avoids small misorientation angles, but a simple geometrical fact that a small 
net misorientation angle demands three independent rotation angles all being small. Indeed, we will see that the 
experimental misorientation distribution vanishes not as 9 2 as would seem natural from equation JJjJ, but as 9. We 
will explain this, and the corresponding cusp in S(n), when we incorporate the effects of cell division, which provides 
a source of new boundaries at zero misorientation angle. 

FigureElshows the misorientations we measure for our Area splitting model. As for the cell sizes, the misorientations 
across cell walls can be averaged weighting them uniformly, by the cell wall perimeter, or by its area. 

We see that the mean misorientation angle in our model grows with e 1 / 2 , as does the experimental misorientations 
across the IDB boundaries. We find similar scaling for Volume splitting, and roughly similar for Uniform splitting 
dynamics; our derivation below suggests that this scaling is generic for our rotational diffusion mechanism. Hence our 
model will not provide an explanation for the 9^ NB ^ ~ e 2 / 3 seen for the GNBs. 

Since the GNBs are in practice distinguished from the IDBs by the number of perpendicular cell walls impinging upon 
them (each GNB has typically a couple of IDBs), we tested whether separating our cell boundaries into previously- 
split GNB analogs and unsplit IDB analogs might lead the misorientations of the former to grow more quickly on 
average with external strain. This did not occur in our model: both previously split and unsplit cell walls scale in 
mean misorientation with the square root of the external strain. 

We can derive a simple differential equation of the time evolution of the distribution of misorientation matrices 
across incidental boundaries 3(n, e). Since new cell boundaries are created at zero angle, the equation will be a 
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FIG. 5: Misorientation Distributions: Area Splitting. The misorientation angles between cell walls, weighting each 
equally (Uniform), or weighting by Perimeter or Area of the cell boundary. LEFT: the mean misorientation angle grows as e 1 ^ 2 
in our model, not only for Area splitting (shown) but also for the other forms of the cell splitting rate. RIGHT: the scaling 
form for the probability distribution is largely independent of the way in which one weights the cells by their size, and agrees 
well with the form p m j a of equation 1141 . 



diffusion equation with a source term: 

fc(n) = 2£V 2 S + Ce&- 1 6(n). (10) 



ck 

where the gradients on the right are with respect to n, and <$(n) is a three-dimensional Dirac delta function (infinite 
at zero, zero elsewhere, integral equal to one). The misorientation matrices diffuse with coefficient 2D because the 
two cells on either side are each diffusing with diffusion coefficient D. 

The first term in equation (|10J) represents the diffusion, or random walk, in rotation space. The second term in 
equation l|l(J|l represents the creation of new cell boundaries which divide the old ones. Because we choose to weight 
our misorientation angle density according to the boundary area, this division does not change the cell boundary 
density 3(n) except at n = 0. 

The new boundary area shows up in our distribution at zero misorientation angle <5(n). To derive the amount of 
new boundary area which is created per unit strain, we use a simple scaling argument. Note that the cell sizes scale as 
L ~ e -1 / 37 . The total boundary area A to t will scale as the number of boundaries 1/L 3 times the area per boundary 
L 2 , hence as e 1 / 37 . The new boundary area needed per unit increment of e thus scales as e~ _1 , giving the prefactor 
for the (5-function in equation i|l(J|) . 

We now specialize to 7 = 2/3 corresponding to Area splitting. If we start with no cell boundaries at e = 0, then 
the solution to equation (|10fl is 

H(n) = / duCu- 1/2 
Jo 

exp(-n 2/ 8L>(e - u))/(8nD{e - u)) 3/2 
= Cexp(-n 2 /8L»e)/(87rL>nVe). (11) 

(The cell boundaries which are formed at deformation u have spread out into a Gaussian of variance D(e—u)) The total 
boundary area A to t — 2Cy / e, as desired. This leads to a prediction for the probability distribution of misorientation 
angles that yields the scaling collapse |QJ 

p(0,e) = OmisWM 

= 6exp(-6 2 /8De)/4:De (12) 

where 



e av = V2irDe (13) 

and 

PmUx) = (ttx/2) exp(— 7re 2 /4) (14) 



8 



(As we will discuss in the appendix „ this happens to be the distribution as derived by Pantleon jlOt lllj in a model 
without cell refinement.) As shown on the right in figure JSJ, this scaling form describes the simulation well: the simple 
scaling argument for the source term of new boundaries (above) captured the behavior of the stochastic simulation. 

One can see from figure that the predicted form is also quite a good description of the experimental data. This 
particular functional form is derived using L ~ e -1 ^ 2 , but the particular exponent is not crucial to the analysis: the 
solution for L ~ e~~ has this general form (with linear slope at 6 = and 9 av (e) ~ yfe) for other 7 as well. 

1 F °° ' ' " " ~" ' : 




FIG. 6: Experiment vs. Theories. Experimental data digitized from Hughes et al.^. The thick curve is p m i a , the scaling 
function of equation Cljl. The thin ClirVG IS /^anisotropic 5 the solution of equation ill ■")! with a length-independent diffusion 
constant with one zero eigenvalue, representing a a highly anisotropic rotational diffusion. 



CONCLUSION: MISSING PIECES 



We believe our simple model is the explanation for the observed scaling seen in cell refinement during plastic flow, 
despite many missing features that are clearly relevant to the experimental behavior. In particular, we would claim 
that any physical system that refines by subdivision, and that undergoes random rotational distortions, is likely to 
enter a scaling regime similar to that seen in our model. In this section, we discuss three of the missing features and 
why they must be relevant for a complete engineering description. 

(1) IDB/GNB Distinction. The experiments show two distinct types of cell boundaries, with different scaling 
exponents and scaling functions. Our model only includes one. While we get a plausible fit to the IDB scaling, 
the GNB misorientation angles grow faster than our rotational diffusion model can reproduce. One must note that 
scaling behavior in other contexts is associated with a single characteristic length scale: our model asymptotes to 
a morphology which is statistically unchanging at large strains except for a single rescaling of length. Because the 
experiments show two length scales, we expect that detailed measurements will show violations of scaling behavior, 
say, when the mean GNB separation length crosses the IDB separation length. Perhaps the fact that the two types 
of boundaries are roughly evolving on perpendicular axes allows for the approximate scaling seen. 

(2) Plastic Deformation. In our simple model, the overall plastic deformation was ignored. In a real material 
under compression the boundaries perpendicular to the strain axis will grow closer to one another linearly in the 
strain even without subdivision! Apparentlv|l2j under the external strain the GNB's slowly rotate towards this axis, 
while the IDB's remain roughly aligned parallel to the axis of compression. The external strain, after this rotation is 
complete, will act to separate the IDB's roughly as the square root of the strain e 1 / 2 . To get the observed refinement 
e -1 / 2 the subdivisions must happen even faster than described by the area splitting law we focused on in the text: 
perimeter splitting, which in the uncle-formed coordinates leads to a length which scales as c -1 would work. This begs 
the question of why the scaling should begin to be observed even before the rotations are complete. 
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Pantleon^^ suggests that this might be an explanation for the observed e -2 / 3 scaling of the separation between 
GNBs. At short times, before the rotation is complete, the lengths vary as e^ 1 ^ 2 , and at long times (with the added 
effects of compression) they will vary faster, perhaps leading to a reasonable fit with the larger exponent. 
(3) Origins of the cell splitting, rotational diffusion. In our model, the cell splitting and rotational diffusion 
are given as part of the dynamics: we do not address what the physical mechanisms are that produce them. 

There are various proposed mechanisms for getting the cell sizes to shrink. Obviously the cell walls cannot just 
move inward: the cell wall velocities would grow linearly in the system size. One could imagine a crinkling of existing 
cell walls (the inverse of the coarsening process seen in spinodal decomposition), or nucleating new cells at junctions 
of existing cell walls: neither picture is compatible with our analysis, and both involve cell wall motion which is 
resisted by pinned sessile dislocation junctions^). Subdivision as we've used it could arise from collisions between 
dislocations as they traverse the cell, although simple calculations suggest that the expected collision rates are too 
small0 E3 and don't scale correctly for our theory 0. Cell splitting due to inhomogeneous stresses induced by 
neighboring cells seem to us the most natural and likely mechanism. A corresponding microscopic picture would 
involve regions in the inhomogeneously stressed cell where the dislocations slow down, or where they are more easily 
pinned by obstacles or other dislocations, leading to the formation of a new cell wall. 

The mechanisms driving the rotations of the crystalline axes of the cells are less well understood. The crystalline 
axes can rotate both directly through the rotation of the material in the cell, and indirectly because of the flux 
of dislocations mediated by the plastic deformation. This latter effect is well studied on larger scales in the field 
of texture evolution, where the plastic deformation of a polycrystal often leads to a gradual alignment of their 
crystalline orientations. Within a single crystalline grain, this texture evolution will on average rotate all the cells 
together: because the relative angles between cells is small, they will largely rotate in the same direction. One 
should note that the traditional explanation for the origin of the GNBs indicates that their misorientation angles 
could well have a overall mean drift in addition to the random diffusion. It is said that GNBs are observed to have 
rotational misorientations that alternate in sign, because they separate regions with differing active slip systems |l3j ; 
the combination of the differing plastic strains and rotations in neighboring pairs of cells can equal the net imposed 
plastic deformation. In this picture, the net rotation angles across the GNB's (equation should have a mean drift 
term in addition to the diffusion term. 
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APPENDIX: CONNECTIONS WITH STOCHASTIC DISLOCATION THEORIES 

The distribution for the misorientation angle for the simple rotational diffusion model, p m is(x), happens to have the 
same form as one derived by Pantleon |lCtllll| without considering cell refinement, and by assuming that the noise in 
the cell orientations was due to random, uncorrelated fluctuations in the dislocation flux. It behooves us, therefore, 
to discuss how Pantleon's work can be interpreted in our context. 

Pantleon's theory, also suggested by Nabarro^j^ and Argon and Haasen0|, is that the stochastic noise in the flux 
of dislocation from either side of the cell boundary leads to randomness in the evolution of the cell boundary angles. 
Each dislocation passing through a cell, say, may shift the top plane of atoms by a distance b with respect to the 
bottom plane, where b is the Burgers vector of the dislocation (roughly the lattice constant). The crystalline axes 
within a cell of characteristic height L will rotate due to one dislocation an amount proportional to b/L. Under a 
strain increment Ae, a cell of characteristic height L must have N = LAe/b dislocations impinging on the side cell 
boundary. A roughly equal and opposite average flux will impinge on the cell boundary from the cell on the other 
side. If the dislocations move independently (which we will argue does not occur), then one expects that there will be 
a net residue after the strain increment of roughly \f~N = LAe/b. Hence the predicted drift in angle after a strain 
increment of Ae is AO = y/Nb/L = \/bAe/L. The diffusion constant is given by D{L) ~ (A9) 2 / Ae ~ b/L. Because 
a single dislocation produces a larger net rotation for smaller cell, the stepsize in the random walk in rotation space 
becomes larger as our cells get smaller. 

Pantleon and Hansen^ consider three cases, where one, two, and three slip systems are activated. The geometry 
of a given cell - the direction and strength of the applied shear with respect to the crystalline axes - will determine 
what types of dislocations are allowed to pass through the cell. If only one slip system is active, the rotation of the 
cell will be confined to a single axis. In our formulation, the diffusion constant in rotation space will be anisotropic: 
in this case it will be a rank one tensor (a 3 x 3 matrix with only one non-zero eigenvalue) . Two slip systems will give 
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a rank 2 tensor, with one zero eigenvalue. Three slip systems will allow the cell to diffuse in any direction, but even 
so the diffusion constant D will in general be anisotropic. 

To make contact with Pantleon, we consider a more general evolution law for the misorientation in area splitting: 



Here we have allowed for an anisotropic, cell-size-dependent diffusion constant by introducing the symmetric tensor 
Dij(L) depending on the current cell size L. The symmetric diffusion tensor Dij will vary with the geometry of 
the particular cell boundary: it can depend both on the relative orientation of the crystal and the strain tensor 
with respect to the plane of the cell boundary. We'll first consider the anisotropy while ignoring the refinement 
(Dij(L) independent of L), in analogy to previous work, and then incorporate the refinement (Dij(L) proportional 
to 1/L ~ e 1 / 2 as argued above). 

Pantleon|3. llCt Hl| in most of his analysis ignores cell refinement. If we make Dij independent of L and set C = in 
equation (|15|l . we get an anisotropic diffusion equation whose solution (assuming a narrow initial distribution of mis- 
orientations) is an anisotropic Gaussian with variances given by the inverse eigenvalues of D, with the experimentally 
observed scaling 8 av ~ e 1 ^ 2 . If we assume D has one non-zero eigenvalue, we get the Gaussian distribution derived 
by Pantleon and Hansen for one active slip system. If we assume D has two equal, non-zero eigenvalues, we get the 
Rayleigh distribution that they find for two perpendicular systems of edge dislocations (which, coincidentally, is the 
same distribution that we found above (|14|l with an isotropic D and a source term). If we assume D is isotropic, we 
get the Maxwell distribution they find for three perpendicular systems of dislocations. 

What happens to the solution of equation l|15fl when we incorporate cell refinement? Pantleon^J notices that 
non-constant cell size must change the scaling of average angle with strain. By setting D{L) ~ e 1 / 2 in accordance 
with reference Q , and changing variables to r = e 3 / 2 in equation i|15|) , we can map it into a problem quite similar to 
equation (|10|l except that the source term is of magnitude proportional to r -2 / 3 . The solution to this equation has a 
shape quite similar to that shown in figure J^J, but with an average ang le 6 av ~ y/T ~ e 3 / 4 . This yields a 3/4 power 
of the strain incompatible with the scaling observed for the incidental cell boundaries. Pantleon is aware of this[lll: 
in an analogous calculation he gets 9 av ~ e ' 72 . 

(If we abandon dislocation noise as the origin of the rotational diffusion, and take the crystalline lattice orientations 
as our basic variables, then there are no microscopic lengths remaining in the problem. The diffusion constant then 
naturally depends only on the local geometry of the cell boundary, and hence is independent of the length scale.) 

The key to the discrepancy is of course that the dislocations motion is not uncorrelated: they must be moving in 
a collective manner. The interaction energy between dislocations is large, and diverges logarithmically with distance, 
reflecting the infinite stiffness of a crystal to gradients in the axes of rotation. In early stages of hardening, one might 
plausibly argue that the dislocations are sufficiently far apart that their interaction forces do not dominate. But in 
the regime studied here, the fact that the dislocations organize into boundaries (avoiding rotational distortions within 
the cells) is a clear signal that it is no longer sensible to treat their evolution independently. If the top half of a cell 
boundary received more dislocations from the right-hand cell than the bottom half, this would produce an enormous 
bending force on the cell. Such an event could only happen for reasonable energy cost if this bending were screened 
by the division of the cell by a new boundary. 



0S(n) 



Dy(i)ViV 3 -3 + Ce 1/2 5(n). 



(15) 



dc 




FIG. 7: Correlated atomic rearrangements associated with simple dislocation motion. 



The stochastic dislocation flux models will be generally applicable whenever the correlations between their motions 
vanish at distances comparable to the cell sizes. Our model is correct in the other limit: our cells only rotate as units. 
Thus our model is appropriate for systems where the dislocation motions are strongly correlated on the scale of the 
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cells - so that their motion can always be described as mediating overall rotations of each cell. Both descriptions are 
only a starting point for a complete theory. 

Consider an analogy: dislocation motion modeled as rearrangements of atoms, versus as the evolution of a continuum 
curve. In the atomic description, one would identify characteristic atomic motions (say, kink diffusion events in the 
case of semiconductors), at rare sites scattered through the crystal. These sites would be strongly correlated (lying 
along the dislocations); their dynamics would produce unusual zipper motions (kinks diffusing along dislocations). 
We expect similar correlations to arise in the dislocation motions mediating cell boundary evolution: the dislocations 
will form correlated dances to keep rotational gradients from entering the cells. On the other hand, treating the 
dislocation motion as the evolution of a continuous curve makes it difficult to incorporate the anisotropic dynamics 
and lattice pinning effects - our model so far has ignored the corresponding crystalline and shear anisotropics in the 
evolution of cell structure. 

We can incorporate some of this asymmetry by hand into our model. Pantleon and Hansen^] point out that the 
individual cell boundaries have rather low symmetries. The diffusion tensor in equation Ijl5|l describes the evolution of 
that subset of cell boundaries with a particular cell boundary orientation and (average) crystal lattice orientation with 
respect to the external shear. There is no reason that for a low-symmetry geometry that the diffusion constant will be 
isotropic, geometry dependent, or material dependent: the general evolution law is i|15[) with Dy independent of L. 
One must solve for the distribution at fixed geometry and then average over geometries to predict the experimental 
distribution (as also discussed in section 5.2 of reference Q). If we assume a strongly anisotropic rotational diffusion 
tensor with one zero eigenvalue and the other two equal (corresponding roughly to Pantleon's analysis with two active 
slip systems) we can solve equation (|15H to find a scaling collapse with scaling function 

P2d\m{x) = (7r 3 cc/64) exp(— ir 3 x 2 /128)K (ir 3 x 2 /128) (16) 

where K$ is the modified Bessel function of the second kind. This scaling function is shown in figure ©. 

It is important also to note that effects that seem clearly related to cell boundary formation have been seen in 
finite-element simulations of single-crystal plasticity by Mika and Dawson et al. 0] . Their system consisted of several 
crystalline grains with differing orientations, subject to an external shear. The inhomogeneous strains within the 
grains led to the formation of subgrain structures very similar to cells. The distribution of cell boundaries in their 
simulation was also found to scale, with strain dependence and functional form similar to that seen for geometrically 
necessary boundaries (that is, scaling with the 2/3 power of strain, and not with the 1/2 power seen for the IDBs). 
Hence Dawson et al. find cells in a simulation totally without dislocations. 



[1] B. Bay, N. Hansen, and D. Kuhlmann-Wilsdorf, Mat. Sci. Eng. A113, 385 (1989). 

[2] D. A. Hughes, D. C. Chrzan, Q. Liu, and N. Hansen, Phys. Rev. Lett. 81 4664 (1998). 

[3] D. A. Hughes, Q. Liu, D. C. Chrzan, and N. Hansen, Acta mater. 45 105 (1997); A. Godfrey and D.A. Hughes, Acta 

mater. 48, 1897 (2000). 
[4] W. Pantleon and N. Hansen, Acta Mater. 49, in press. 

[5] J. P. Sethna and E. Demler, cond-mat/0104552 version 2 (early draft of current manuscript). 

[6] D. P. Mika and P. R. Dawson, Acta Materialia 47, 1355 (1999). N. R. Barton and P. R. Dawson, "On the spatial 

arrangement of lattice orientations in hot rolled multiphase titanium", to be published. 
[7] E. D. McGrady and R. M. Ziff, Phys. Rev. Lett. 58, 892 (1987). 
[8] Z. Cheng and S. Redner, Phys. Rev. Lett. 60, 2450 (1988). 

[9] P. Hahner, K. Bay, and M. Zaiser, Phys. Rev. Lett. 81, 2470 (1998), R. Thomson and L. E. Levine, Phys. Rev. Lett. 81, 
3884 (1998), I. Groma and B. Bako, Phys. Rev. Lett. 84, 1487 (2000). 
[10] W. Pantleon, Scripts Materialia 35, 511 (1996). 
[11] W. Pantleon, Acta Mater. 46, 451 (1998). 
[12] W. Pantleon, private communication. 
[13] D. Hughes, private communication. 

[14] F. Prinz and A. S. Argon, Physica status solidi 57, 741, (1980). They discuss not collisions between dislocations on different 
slip systems, but mutual trapping of dislocations with opposite Burgers vector on one slip system: our analysis still applies. 
One should note that Argon now believes that cell boundaries may move during plastic deformation through avalanches 
of failures of these sessile dislocation junctions. 

[15] F. R. N. Nabarro, Scripta Metallurgies et Materialia 30, 1085 (1994). 

[16] A. S. Argon and P. Haasen, Acta Metall. Mater. 41, 3289, (1993) equation 21. 

[17] Write V = J p vo \{V)dV, use equation [5] to compute dV/da, change variables to y — V/V, and substitute in the scaling 
form equation |21 



